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ABSTRACT 

We present the results of H- and K-band VLT/SINFONI integral field spectroscopy 
of the ULIRG IRAS 19254—7245 (The Super-antennae), an interacting double galaxy 
system containing an embedded AGN. Deep i^-band spectroscopy reveals Paa arising 
in a warped disc with position angle of 330° and an inclination i = 40 — 55°. The 
kinemetric parameters derived for II2 are similar to Paa. Two high- ionization emission 
lines, [Sivi] and [Alix], are detected and we identify as [Nill] the line observed at 
1.94 /im. Diluting non-stellar continuum, which was previously detected, has decayed, 
and the JY-band continuum emission is consistent with pure stellar emission. Based 
on II2 emission line ratios it is likely that at the central 1-kpc region H2 is excited by 
UV fluorescence in dense clouds while shock excitation is dominant further out. This 
scenario is supported by very low Paa to II2 line ratio detected outside the nuclear 
region and non-thermal ortho/para ratios (~2.0 - 2.5) close to the nucleus. 

Key words: galaxies: starburst - galaxies: Seyfert - galaxies: individual (The Super- 
antennae) - infrared: galaxies 



1 INTRODUCTION 

The Super-antennae (IRAS 19254-7245) at a redshift of 
z ~ 0.062 (distance of ^ 260 Mpc assuming Hq — 71 km s~^ 
Mpc"^ Spergel et al. 2003) is an interacting double galaxy 
system where the two components are separated by ^8^^('-^2 
kpc). Visually the most notable feature of the system are 
the thin antennae extending up to 5^ ('^350 kpc), far longer 
than similar feature of the Antennae galaxy (^100 kpc). 

The starburst activity powering the total infrared lumi- 
nosity of 10^^-^^ L/i^/L© (Due, Mirabel & Maza 1997) in the 
system has been triggered by a 3 : 1 mass ratio encounter 
(Dasyra et al. 2006). The northern nucleus is less luminous 
than the southern nucleus in the infrared, and is believed to 
be in a post-starburst stage (e.g. Berta et al. 2003), while 
the southern galaxy is a starburst with an embedded AGN. 
Berta et al. (2003) found the spectrophotometric character- 
istics of the southern nuclei to be mostly consistent with two 
main stellar populations: an intermediate aged (IGy) star- 
burst, which represents '^35 percent of the luminous mass 
of the galaxy and an old 12 Gy population representing the 
remaining ^-^65 percent of the total luminous mass. In addi- 



Based on observations collected at the European Southern Ob- 
servatory, Paranal, Chile (60A-9041(A)). 



tion, the southern nucleus has an intense ongoing starburst, 
which contributes little to the total luminous mass. 

The evidence for the presence of an AGN comes from 
the optical emission line ratios (Mirabel, Lutz & Maza 1991), 
IR coronal-emission lines (Vanzi et al. 2002), observations 
of a hard, flat X-ray component above 2 keV (Pappa, Geor- 
gantopoulos & Stewart 2000), and strong Fe Ka emission 
lines (Braito et al. 2003) interpreted as a Compton-thick X- 
ray source. Like in other LIRGs/ULIRGs, a long-standing 
question has been whether the spectral energy distribution 
of the southern galaxy is dominated by the starburst or 
the AGN. The mid-IR spectroscopic diagnostics are inter- 
mediate between those for starbursts and AGN (Genzel et 
al. 1998). Recently, Charmandaris et al. (2002) and Risal- 
iti et al. (2003; 2006) claimed the AGN is the main energy 
source in the southern galaxy based on mid-IR ISO and deep 
L-band ground-based spectroscopy, respectively. 

Based on the integral field spectroscopic data obtained 
with VLT/SINFONI on the southern galaxy, we discuss the 
nuclear spectrum, excitation of molecular H2 emission lines 
and the morphology and kinemetric properties of ionized 
and molecular gas. 

In addition to the high- excitation emission lines requir- 
ing the presence of either strong ionising continuum or fast 
shocks with Vs — 200 km s~^ or greater (Contini & Viegas 
2001), other AGN indicators in near-IR wavelength range 
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covered in this work are the variabihty of the emission hnes 
and non-stellar continuum. Considering the distance and 
brightness of the Super-antennae, supernovae and their sub- 
sequent evolution are unlikely to produce detectable spec- 
troscopic variability at the nucleus. 

The dominant excitation mechanism for H2 has been 
under debate for years, the two main candidates being UV 
fluorescence (e.g. Black & van Dishoeck 1987; Sternberg & 
Dalgarno 1989) and shocks (e.g. Hollenbach & McKee 1989). 
The interpretation of H2 line ratios is complicated due to 
the strong dependence of line ratios on density; the energy 
levels are thermalised through inter- molecule collisions and 
the line ratios in dense clouds illuminated by intense UV 
radiation are very similar to those from purely thermal ex- 
citation, either through shocks or X-rays. Furthermore the 
H2 line ratios are affected by the spin degeneracy of the lev- 
els. Radiative decay between ortho (odd J) and para (even 
J) is not possible due to different spin and in the ground 
electric state only A J = 0, ±2 transitions are possible. For 
purely thermal excitation the ratio between the ortho and 
para levels is 3.0, but is lower for fluorescence in tenuous gas 
(-1.9 using the 1-0 S(l), S(2) and S(3) transitions; Black 
& van Dishoeck 1987). 

As the Super-antennae is an interacting system, the 
kinematics is expected to be more complex than in the case 
of non-interacting galaxies. In the literature multiple gaus- 
sian have been used to approximate optical emission lines 
(Vanzi et al. 2002; Colina, Lipari & Macchetto 1991). In- 
terpretating the results of multi-gaussian fits in a physically 
meaningful manner is however non-trivial, unless the com- 
ponents are clearly separated from each other. A more qual- 
itative approach for tracing multiple velocity components is 
using the Gauss-Hermite series to describe the line profile 
(van der Marel & Franx 1993; Gerhard 1993). Strong kine- 
matical components can also be revealed by fitting a warped 
disc model to the observed velocities. 

The paper is organised as follows. In Section 2, the ob- 
servations and data reduction are described. In Section 3 
the results are presented and in Section 4 we summarise the 
conclusions of the work. 



2 OBSERVATIONS AND DATA REDUCTION 

The southern galaxy of the Super-antennae was observed 
during SINFONI (Spectrograph for INtegral Field Observa- 
tions in the Near Infrared; Eisenhauer et al. 2003; Bonnet 
et al. 2004) Science Verification in Paranal, Chile, in August 
2004 both in H {Rc^ 2900) and K bands {R 4500) with 
a 2048^ pixel engineering grade detector. This engineering 
grade detector had two significant defects: a "glowcenter" 
at 2.38 /xm and a large dead region at >2.4/im. Unfortu- 
nately these artifacts compromised the data quality at the 
nucleus. The 250 mas pixel scale was used during the obser- 
vations without Adaptive Optics providing a field of view 
(f.o.v) of 8^^x8^^ The total integration time was 14000 s in 
K and 1800 s in if- band. The average l/-band seeing mea- 
sured by the Differential Image Motion Monitor (DIMM) 
during the observations was 0^.^94 and 0^.^77 during the K- 
band and if-band observations, respectively, and the con- 
ditions were non-photometric. The observations were done 



following the classical OSSO sequence of alternating object 
and sky exposures. 

The data reduction was done with the spred software 
package (Schreiber at al. 2004) developed at the Max Planck 
Inst it ut fiir Extraterrestrische Physik. The data were bad 
pixel and distortion corrected and sky subtracted using the 
nearest sky exposure. Instrumental flexure was tracked by 
cross-correlating the night sky OH emission lines before sky 
subtraction, and wavelength calibration was done by fit- 
ting a quadratic polynomial to Ne arc lamp frames for the 
K-hand and Xe+Ar frames the if-band. Due to unstable 
conditions during the observations, the night sky emission 
lines and thermal background at the red end of the K band 
were not perfectly removed. This residual background was 
removed by calculating a median value along the edges of 
the cubes at each spectral plane and subtracting it. While 
this subtraction unavoidably removes faint continuum emis- 
sion from the galaxy, the emission lines are unaffected as not 
even Paa can be detected this far out and is not present in 
the subtracted spectra. 

The relative wavelength shifts caused by the instru- 
mental flexure between different exposures were determined 
by cross-correlation after collapsing the reconstructed data 
cubes. The cubes were combined and divided by the telluric 
standard spectra (Hip099481; B9.5V class star) in order to 
remove atmospheric signatures. Finally, the flux calibration 
was done by comparison with broad-band images (Vanzi et 
al. 2002). We estimate the absolute flux calibration is accu- 
rate to within 5-10 percent. 

In order to enhance the signal-to-noise ratio of the im- 
ages the data were smoothed with a 3 pixel x 3 pixel boxcar 
(0^.^38 X 0^.^38), which is much smaller than the resolution of 
the data. Emission line measurements were done automati- 
cally with scripts written by us for IRAfQ- If a detected line 
was significantly narrower than the instrumental profile, it 
was subtracted from the spectra and the fit was repeated. Fi- 
nally, the detections were inspected and all obviously wrong 
fits due to cosmic ray residuals or noise were rejected. In 
the case of Paa, which has multiple components that are 
difficult to fit automatically, the fit was carried over inter- 
actively within the central region. In all the images shown 
in later sections the lowest level displayed is at 2.5cr unless 
otherwise stated. 



3 RESULTS 

A broadband F804W image of the Super-antennae obtained 
from the HST data archive is presented in Fig. [1] The SIN- 
FONI field was centred on the southern nucleus and the field 
of view contains part of the bridge detected in broadband 
imaging connecting the two nuclei. 

3.1 Nuclear spectrum 

The nuclear spectrum of the Super-antennae is displayed 
in Fig. [2] In addition to Paa fortuitously located at an at- 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 2. The summed spectrum of the central 5x5 spectra (0(^6x0^/6) in units of lO"-*^^ ergs cm"-*^ s"-*^ /im"-*^ with hne identifications 
as indicated {upper panel). In the lower panel the normahsed residuals after subtracting the strongest Hi (Paa, Br7) and H2 lines (1-0 
S(l), S(2), S(3)) and smoothing with residual with five pixel boxcar are shown. The x-axis represents the observed wavelength. 
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Figure 1. HST F804W filter image of the Super-antennae. The 
SINFONI field of view is indicated with a box. The antennae of 
the galaxy do not fit within the area depicted in this figure. East 
is left and north is up in this and any subsequent images. 



mospherically relatively clean region at ~2 /xm, the i^-band 
spectrum shows many prominent molecular H2 lines, hydro- 
gen recombination lines Br7 and BrJ and Hex 2.058 /xm. A 
careful study reveals also two coronal lines - [Si Vl] 1.962 /im 
and [Mix] 2.04 /xm (redshifted to 2.0833 /xm and 2.166 /xm, 
respectively). We have compared the observed fluxes to the 



ones reported by Vanzi et al. (2002) in identical apertures. 
While the H2 line ratios are in good agreement with Vanzi 
et al. (2002), the hydrogen recombination lines Br7 and Br(5 
are much weaker relative to H2 lines. As the H2 lines are 
unlikely to change in time scales of a few years, this vari- 
ability in line ratios strongly suggests large fraction of the 
hydrogen recombination line luminosity in near-IR is due to 
AGN illuminated gas in NLR. The H2 emission line fluxes 
are higher by ~40 percent than those reported by Vanzi et 
al. However, as discussed above, we were forced to use the 
broad-band images of Vanzi et al. for flux calibration, and 
it is possible that this difference can at least partly be at- 
tributed to changes in continuum brightness. This interpre- 
tation is likely, as the if-band CO absorption lines are less 
diluted by non-stellar continuum (Sect. 13. 3p than detected 
by Vanzi et al. (2002). 

A clear indication of an AGN in the Super-antennae is 
the presence of coronal emission lines: [Si Vl] (ionisation en- 
ergy 167 eV) previously detected by Vanzi et al. (2002) and 
a faint [Alix] emission line (285 eV) detected at 12cr level. 
While [Al ix] has been detected only in few AGN (Maiolino 
et al. 1999; Reunanen et al. 2003), the line is one of the 
strongest NIR coronal lines when the ionising parameter 
(the ratio of hydrogen ionising photon to total hydrogen 
densities) is high enough. Detection of highly-ionised emis- 
sion lines is not typical in ULIRGs; Murphy et al. (2001) 
found [Si vi] in only two galaxies out of 33 ULIRGs in their 
survey. This simply reflects the fact that stars are unable 
to photo-ionise the coronal emission lines, and the shocks 
in starforming regions are not likely excitation mechanism 
(e.g. Marconi et al. 1994; Rodrfgues-Ardila et al. 2002). 

An intriguing similarity between our nuclear spectrum 
and that presented in Vanzi et al. (2002) is the presence 
of a faint emission line at 2.058 /xm, which they attributed 
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Figure 3. The Pao; spectrum summed over the central 0(^6x0(^6. 
In the upper panel the fits with one (long- dashed), two (dashed) 
and three Gaussians (dotted) are indicated as are the H2 emission 
lines. In the lower panel the residuals of the fits are shown. The 
shape of the H2 lines is taken from the isolated 1-0 S(l) line. The 
inset shows the comparison between the three component model 
of Ho; (Vanzi et al. 2002) and Paa. 



Table 1. Observed emission lines in the nucleus. Fluxes are 
given in units of lO"-*^^ erg s"-*^ cm~^. [Fell] 1.64 /xm is blended 
with Brackett 12-4 1.64117 /im, [Fell] 1.68 /xm with Brackett 11-4 
1.68112 /im and BrS with H2 2-1 S(5) 1.94487 /im. 



Line 



Ao 



A. 



Flux 



FWHM 

km s"-*^ 



[Fe 11] 


15338.9 


16314.3 


±1.4 


1.09 ± 0.06 


1890 ± 80 


[Fe 11] 


16439.9 


17457.3 


±0.4 


5.10 ± 0.07 


1110 ± 10 


[Fe 11] 


16773.3 


17826.4 


±3.0 


1.24 ± 0.10 


1950 ± 130 


1-0 S(5) 


18358.0 


19472.2 


±3.1 


2.64 ± 0.4 


1270 ± 160 


Paa 


18756.3 


19909.8 


±0.2 


13.30 ± 0.04 


515 ± 2 


1-0 S(4) 


18919.7 


20090.8 


±1.8 


0.52 ± 0.05 


830 ± 80 


[Ni 11] 


19393.0 


20586.9 


±2.2 


0.16 ± 0.04 


550 ± 110 


Bid 


19451.0 


20652.1 


±1.1 


0.37 ± 0.03 


600 ± 40 


1-0 S(3) 


19575.6 


20798.2 


±0.3 


3.34 ± 0.03 


920 ± 10 


[Si VI] 


19634.1 


20899.3 


±0.9 


0.72 ± 0.02 


900 ± 30 


1-0 S(2) 


20337.6 


21609.2 


±0.3 


1.20 ± 0.02 


1140 ± 20 


[Alix] 


20449.9 


21695.6 


±1.3 


0.12 ± 0.01 


450 ± 50 


He I 


20586.9 


21863.7 


±0.6 


0.25 ± 0.01 


550 ± 30 


2-1 S(3) 


20735.1 


22020.4 


±1.1 


0.13 ± 0.02 


500 ± 50 


1-0 S(l) 


21218.3 


22536.5 


±1.1 


2.99 ± 0.02 


953 ± 4 


Br7 


21661.3 


23002.3 


±0.6 


0.66 ± 0.02 


540 ± 20 


1-0 S(0) 


22233.0 


23607.7 


±1.0 


0.40 ± 0.03 


580 ± 40 



to night sky OH emission lines. However, the same feature 
is also present in our spectra which is otherwise free of tel- 
luric night sky signatures. The presence of several, relatively 
broad emission lines (2-1 S(5) 1.94487 /xm, Br^ 1.9451 /xm, 
1-0 S(3) 1.95756 /xm, [Si vi] 1.96341 /xm) makes the study of 
the region complicated. In order to reduce the number of 
free parameters in the fit, we have assumed the line profiles 
of 1-0 S(3) and BrS are identical to those of the isolated 1-0 



S(l) and Paa Unes, respectively. Furthermore, we have used 
the 1-0 S(l) line profile for the 2-1 S(5) line fitting, keep- 
ing the position of the line fixed relative to 1-0 S(3) line. 
The addition of the 2-1 S(5) line, which is strongly blended 
with Br^, is necessary to make the residuals around Br^ 
smaller. In typical starforming galaxies with 1-0 S(l)/Br7 
-1 (e.g. Goldader et al. 1997), the ratio 2-1 S(5) to Br^ 
is ^0.06 assuming thermal excitation with Tex — 2000 K. 
However, due to the rather extreme 1-0 S(l)/Br7 ratio of 
'^5 observed in the Super-antennae, the contribution of 2- 
1 S(5) is ^30 percent of Br^ flux. The inclusion of 2-1 
S(5) does not affect the position of the unidentified line. 
The fitting gives the wavelength for the unidentified line 
2.0584±0.0001 /im (Ao ^ 1.937 /im). The best candidate in 
the NIST Atomic Spectra Database is [Niii] 19393A corre- 
sponding to the 3d^ 4s a^F-3d^ 4s a^F transition. Optical 
emission lines (AA 7380 and 7414) originating from the same 
upper energy level have been detected in objects of various 
types, including nebulae (e.g. Osterbrock, Tran & Veilleux 
1992), supernova remnants (e.g. Dennefeld 1986) and AGN 
(e.g. Halpern & Oke 1986). As the 1.93 /xm line should also 
be relatively strong (based on the relative A- values), the 
identification of this line as [Ni il] seems secure. The optical 
lines have, however, not been detected in the Super-antennae 
(Berta et al. 2003). 



3.2 Emission line profiles 

The most notable feature of the nuclear spectrum is of course 
Paa. However, Paa is badly fitted with a single Gaussian 
(Fig. [3]). Similarly, multiple components (three or more; 
Vanzi et al. 2002; Colina, Lipari & Macchetto 1991) have 
been used to approximate the optical emission lines. The 
wavelength region (^1.99 /xm) around Paa appears to be 
completely free of atmospheric (telluric) residuals - redward 
of Paa the strength, shape and position of 1-0 S(4) at 1.89 
/im is in excellent agreement with the measurements of the 
stronger and isolated 1-0 S(l) Une. Therefore, the contin- 
uum level around Pa line is very well determined. Paa is 
asymmetric, with the blue wing being stronger, and is sim- 
ilar to the Ha profile observed by Vanzi et al. While a two- 
component fit produces acceptable results in the sense that 
the residuals (red asymmetry) can be readily explained by 
obscured broad Une region, the similarity with the Ha sug- 
gests this is not the case. However, we note that the FWHM 
of the broader component in the two-component fit is simi- 
lar to possible detection of broad Ha in polarised emission 
by Pernechele et al. (2003). 

When compared to Ha, the center- most component in 
Paa is much stronger and represents ^50 percent of the 
total Paa emission, while the narrow (a = 185 km s""*^) 
component of Ha provides only 6.5 percent of the total line 
emission (Vanzi et al 2002). Qualitatively the weakness of 
coronal emission lines relative to the hydrogen recombina- 
tion lines when compared with other Sy2 galaxies suggest 
the central component arises from a starburst and not from 
the NLR. 

All the three main emission lines (Paa, H2 1-0 S(l) 
and [Fell]) have different line profiles. H2 has a strong red 
component, while [Fell] has a weak asymmetry to the red 
and is broad. Since [Fell] has critical density of the order 



of iVe 



10^ 



10^ cm""* (Blietz et al. 1994), the BLR can- 
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not contribute significantly to the flux. Thus the shape of 
[Fell] also supports multiple velocity components instead 
of a genuine broad component. Careful deblending also sug- 
gests that [Si vi] may be double-peaked, with a separation of 
-750 km s"^ {z = 0.06155, 0.06419) between the two peaks, 
indicating it may originate from an unresolved disk. [Alix] 
is too weak to determine whether the profile is similar to 
[Si VI]. 

3.3 Nuclear absorption lines 

In addition to many prominent emission lines, the nuclear 
spectrum also contains several absorption lines, allowing the 
application of the penalised pixel fitting method (Cappel- 
lari & Emsellem 2004) to determine the line of sight ve- 
locity distribution. In this process, the line of sight veloc- 
ity distribution is described by a Gauss-Hermite series (van 
der Marel & Franx 1993) and the solution is based on the 
maximum penalised likelihood formalism (see Cappellari & 
Emsellem 2004 and references therein for details). The en- 
gineering grade detector used in SINFONI during the time 
of the observations had a large dead region at the red end 
of the spectrum, which compromised the data quality in the 
i^-band beyond >2.4 jim.. Therefore, only the if-band was 
used. The fitting was done using only the low-order moments 
{v^ cr), as including the higher orders hs and did not pro- 
vide a significant improvement in the quality of the fit, most 
likely due to relatively low S/N and beam-smearing. Late 
type stars (K3V, MOV, M5III) obtained with SINFONI for 
another observing program were used as templates. 

We derive a — 210 ± 20 km s~^, in agreement with a — 
175 ± 24 km of Dasyra et al. (2006) and 188 ± 10 km 
of Tacconi et al. (2002). The H-band continuum is consistent 
with late type stellar template, with no need to introduce 
a non- stellar component to the fit. This is in contrast to 
the findings of Vanzi et al. (2002), according to whom the 
non-stellar component contributed 30 percent of the if-band 
emission based on a similar fitting on the if-band and an 
even higher fraction based on broad-band colours. 

The disappearance of the non-stellar component is not 
surprising as our observations and those of Vanzi et al. were 
separated by 5 years. The most likely origin of the dilut- 
ing continuum is hot dust heated close to the sublimation 
temperature, and therefore follows closely changes in the 
heating flux. Alternatively, the diluting continuum may be 
due to synchrotron, as has been suggested for some of the 
nearby AGN in near-IR (Prieto et al., in preparation). In 
both cases the implication is that, at least in near-IR, the 
non- stellar continuum is solely due to the AGN. As discussed 
above, the Super-antennae may be AGN dominated even in 
mid-IR (Charmandaris et al. 2002; Risaliti et al. 2003; 2006). 

3.4 Emission line morphologies 

The emission line fluxes relative to the nuclear peak flux, 
the full width at half-maximum (FWHM) of the Gaussian 
fitted to the line profile and the velocity field for Paa 1.87 
/im, [Fell] 1.64 /xm and the H2 line 1-0 S(l) 2.12 /xm are 
displayed in Fig. 2] For comparison, the HST 0.8 /im im- 
age mildly smoothed to 0^.^3 resolution is shown in contours. 
The alignment between the HST image and SINFONI data 



was based on the position of the nuclear continuum source. 
While the southern nucleus is red, and the optical nucleus 
may therefore not coincide with the near-IR nucleus, no ev- 
idence for this is found based on comparison between large- 
field near-IR and HST images, at least within the seeing- 
limited resolution of data discussed in this paper. In case of 
Paa, the continuum wavelength range in automatic line fit- 
ting has been selected to exclude contribution from "broad" 
component; thus the fitting procedure (and line ratios) is 
comparable with the other lines where the Paa components 
cannot be detected. 

As the SINFONI field of view does not contain any 
stars, a direct estimate for the achieved spatial resolution 
cannot be derived based on the data itself. Furthermore, the 
DIMM sensor measures optical seeing close to zenith while 
the Super- antennae was observed at an airmass of '-^1.5. 
The nucleus is unresolved in the iC-band image (Vanzi et 
al.), where the spatial resolution ^^l^.^O (1.3 kpc) can be de- 
rived directly from the numerous stars present in the large 
f.o.v. After collapsing the spectral cubes in emission line- 
free wavelengths, we derive the FWHM size for the nuclear 
continuum source 0^9 both in H and similar to that mea- 
sured from pure emission line maps for the main emission 
lines (Paa, [Fell], H2 lines). While these values are larger 
than reported by the DIMM in the optical, especially dur- 
ing if-band observations, the derived nuclear sizes are in- 
consistent with the unresolved nucleus in broadband images. 
More likely, the DIMM sensor has underestimated the opti- 
cal seeing and the nucleus is unresolved in both the collapsed 
spectral continuum images and emission lines. 

While the nuclear source itself is unresolved, a fainter 
extended component can be detected in most emission lines. 
The most extended emission line is Paa 1.8756 /xm, which is 
dominated by the nucleus with spiral structures extending 
to south-east and north-west. These Paa features coincide 
with the dusty spiral present in HST archive images (Fig.[T]). 
Similar to Ha emission (Mihos & Bothun 1998), both Paa 
and molecular H2 lines show enhanced emission northwest of 
the nucleus. The southern spiral consisting of several clumps 
(starforming regions) is also faintly traced in H2. Further- 
more, a separate starforming clump is seen only in Paa 
north of the nucleus coinciding with the bridge detected in 
continuum (Fig. [1]) connecting the southern galaxy to the 
northern. 

3.5 Molecular gas 

Davies et al. (2003) studied the molecular gas in 7 ULIRGS 
and found that while the 1-0 transitions appear to be ther- 
malised, the same is not true for 2-1 and 3-2 transitions. 
They concluded the H2 excitation is predominantly UV flu- 
orescence, but the lower levels are thermalised by inter- 
molecule collisions in dense clouds. This produces different 
excitation temperatures, ~1300 K based onv — 1 transitions 
and '^5000-6000 K based on v — 2 and v — 3 transitions. 
In order to probe the fainter transitions, we have corrected 
the spectra for rotation with velocities derived from the 1- 
S(l) line. This "derotation" increases S/N of the spectra 
by decreasing the width of the spectra features and decor- 
relates possible telluric residuals; this procedure is identical 
to what was used in Davies et al. (2005). However, due to a 
lack of CO-bandheads at iT-band, the template fitting is so 





Figure 4. From left to right, the flux normahsed to the peak nuclear flux, Gaussian FWHM and the velocity fleld for Pao; (top row; 
excluding the "broad" component), [Fell] (middle row) and H2 1-0 S(l) (bottom row). The optical morphology from HST 0.8 /im image 
is shown in contours. 
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Table 2. Comparison between the H2 emission line ratios in 
ULIRGs, AGN and the Super-antennae; the Br7 for the extended 
case is estimated from Pao; 



Hg ortho/para 



Line 


ULIRG 


AGN 


Super-antennae 








Nucleus 


Extended 


1-0 S(3) 






0.595 


0.739 


1-0 S(2) 


0.337±0.018 


0.379±0.068 


0.496 


0.327 


2-1 S(3) 


0.146±0.022 


0.217±0.132 


0.045 


0.045 


1-0 S(l) 


1.000 


1.000 


1.000 


1.000 


2-1 S(2) 


0.074±0.007 


0.118±0.012 


0.039 


< 0.06 


1-0 S(0) 


0.271±0.035 


0.327±0.047 


0.308 


0.350 


2-1 S(l) 


0.142±0.031 


0.125±0.038 






Br7 


0.6-1.33 




0.256 


0.164 



badly constrained, that we are unable to remove the stellar 
contribution. 

The comparison between the average line ratios pre- 
sented by Davies et al. (2003) for ULIRGS, Davies et 
al. (2006) for AGN and the Super-antennae is presented 
in Table (2] For the Super-antennae, the line ratios both 
for the nuclear region (r <0^/6) and extended emission 
(IV25 < r < 2V5) are given. When compared with the 
ULIRG sample, the 2-1 transitions appear weaker both at 
the nucleus and in the extended annulus, and are not over- 
populated relative to the 1-0 transitions. Br7/1-0S (1) is 
much lower at the nucleus, and in the extended annulus Br7 
cannot be even detected. Assuming the Case B Paa/Br7 ra- 
tio of 12.1, we infer a Br7/l-0 S(l) ratio of -0.16 there. The 
decrease in Paa/1-0 S(3) ratio is also clearly visible in Fig. 
[6l which interestingly also reveals a small ring-like structure 
with diameter —1^5. 

We have estimated the excitation temperature from 1-0 
S(3)/l-0 S(l), and ortho/para ratio from 1-0 S(2)/l-0 S(l) 
since these are the strongest H2 emission lines in the spectra. 
The ortho/para ratio can easily be shown to be practically 
independent of extinction: 



fo 



Iv 



S{2)/S{1) 



-10 



-0.004 Aj^ 



(1) 



where S{x) denotes the observed intensities for 1-0 S(x) 
emission lines and Ak is the K-hand extinction assuming 
Ax (X A~^"^^ and the A-coefficients from Wolniewicz, Sim- 
botin & Dalgarno (1998). The ortho/para map is shown in 
Fig. \5\ showing a relatively constant ratio. By summing the 
"derotated" spectra in annuli, we can follow the ortho/para 
ratio changes further out than calculating it in individual 
spectra. At the nucleus (r < 0^6) the ratio is 2.4zb0.1, and 
remains similar in the annulus 0.6 < r < 1.25^^ Further 
out, the ratios are consistent with thermal excitation of H2 
(ortho/para ^^3). As the minimum values detected in or- 
tho/para ratio coincide with the ring- like structure visible 
in Paa/1-0 S(3) ratio, it is likely that in the nuclear re- 
gion H2 is excited by UV fluorescence in dense clouds. A 
contribution from X-ray excitation may be needed to de- 
populate the 2-1 states (e.g. Gredel & Dalgarno 1995), as 
Davies et al. (2003) observed larger ratios in their ULIRG 
sample. Similar scenario has been suggested for NGC6240 
(Draine & Woods 1990), though higher- resolution data is 
clearly needed to properly resolve the possible starbursting 
ring and separate the mechanisms. 
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2 
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Figure 5. Ortho-Para ratio of H2 determined from 1-0 lines S(l), 
S(2), S(3) 



Paa / Hg 1-0 S(3) 




2 

R.A. [arcsec] 



1.5 2 2.5 3 3.5 



Figure 6. The Paa/H2 1-0 S(3) ratio 



For the extended H2 gas the situation is more clear, as 
both the ortho-para ratio and low Paa/H2 line ratio suggest 
thermal excitation through shocks. While in principle X-rays 
can also produce low Paa/H2 ratios (Gredel & Dalgarno 
1995), such a scenario seems unlikely, as the likely origin for 
the extended X-ray emission is supernovae; however, [Fell] 
is strongly concentrated at the nuclear region. 

The average column density of molecular hydrogen 
within the 5x5 pixel region is 1.8 x lO^^cm"^, correspond- 
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ing to the total mass of hot excited molecular hydrogen of 
1.8 X 10^ M0. The peak (seeing dependent) column density is 
a factor of 1.3 bigger than the average value. The total mass 
of the available molecular material as determined from CO 
radio emission is several orders of magnitudes larger, 3x 10^° 
M0 (Mirabel et al. 1990); however, the radio beamsize (44") 
includes both the nuclei in the system. 

3.6 Kinematics 

Low resolution Ha velocity maps for the Super-antennae 
were presented by Mihos & Bothun (1998), whose Fabry- 
Perot images also covered the northern galaxy. The most 
notable difference between our maps and theirs is the plume 
detected in Ha in East of the nucleus, which is not very 
extended in Paa. This may be simply due to lower signal to 
noise ratio at the outer edges of the mosaiced spectra. The 
faint streaming motion with Av —130 km~^ appears to 
be present in Paa in South of the nucleus extending up to 
4" (Fig.lU). 

We have applied the Kinemetry IDL routine described 
by Krajnovic et al. (2006) to the data. This method is a 
generalisation of surface photometry performed in isophotal 
annuli of galaxies applied to higher order moments (v, a, /i3, 
/i4...)- The routine can either use a constant systemic veloc- 
ity or allow it to change as a function of radius; for isolated, 
non-interacting galaxies the systemic velocity (Aq follow- 
ing the terminology in Krajnovic et al. 2006) is expected to 
be nearly constant. We have tried both methods with Paa 
in the Super-antennae, but constant Aq produces a poorer 
match with the observed velocity field. The parameters of 
the warped disk model with Aq kept as a free parameter are 
shown in Fig. [71 For the four innermost radii, the fitting is 
done inside the PSF and therefore the position angle and 
ellipticity have very large errors. This simply represents the 
fact that while the fitting process is numerically stable even 
inside the HWHM radius of the PSF, the shape parameters 
derived there have no real physical meaning. 

The fitting has revealed the presence of at least two 
separate kinematic components within the central r — 4"'m. 
the Super-antennae, possibly three. In the inner parts of the 
galaxy the velocity field is well modelled by a warped disk 
model, and the residuals between the observed velocities and 
model are generally within ^25 km s~^ . Northeast and west 
of the nucleus, the residuals are much larger, up to 150 km 
s~^. While the shape parameters [q and position angle) of 
the outermost ellipses does depend on how large the area 
the fitting is carried out over or whether the beginning of 
the spirals arms is included or not, the closeness of the edge 
of usable data does not affect these residuals significantly. 
Especially northeast of the nucleus, the area of large residu- 
als is also clearly visible in the observed velocity field (Fig. 

a. 

The ellipticity of the fitted rings, related to the fiat- 
tening parameter of Krajnovic et al. (2006) by ^ = 1 — e, 
gradually increases from 0.24 at r = 0'.'68 to 0.45 at r > 2" , 
while the position angles of the ellipses are almost constant, 
327-331° . Aq represents the systemic velocity of the fitted el- 
lipses and is relatively stable between 0^.^8 - 2" ^ '-^18530 km 
s~^. Further out, Aq decreases and the ratio increases. 
As the presence of multiple components is detected by the 
rise of the coefficient in parametrisation of Krajnovic et 
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Figure 7. Kinematic model parameters for the Pao; emission 
{circles with error bars) and H2 1-0 S(l) (open squares). From 
top to bottom: the position angle of the ellipses, the flattening 
coeflicient q (related to e by g = 1 — e), the systemic velocity Aq, 
rotational velocity ki and non-rotational velocity component k^. 



al. (2006), this indicates that the outer parts of the galaxy do 
not follow the same disk-like kinematics as the inner parts. 
This outer component is already clearly visible in the ob- 
served velocity field (Fig. |4|. Further support for the sepa- 
rate outer component comes from the Gauss-Hermite expan- 
sion of the Paa profile (Fig. [5]), where the higher order mo- 
ments hs and /i4 are directly observable quantities without 
any modelling assumptions, /is and /i4 in the inner regions 
closely follow the expectations for a warped disk located at 
position angle of ~330°. Northeast and west of the nucleus, 
both hs and /i4 show significant deviations from disk-like 
rotation; these deviations coincide with the the regions with 
Vobs — Vmodei < "25 km s""*^ of the model fitting. 

Closer still to the nucleus, with r < 0^7, the fitted Aq 
values are higher than in the surrounding region. The na- 
ture of this jump is not clear based on the current, seeing 
limited data. The redshifted nuclear velocities appear to be 
connected with the positive residuals present in Vobs—Vmodei 
map at position angle of ~70°(Fig. [8]). 

We have repeated an identical fitting procedure with 
1-0 S(l) as we did with Paa. Overall, the quality of the 
kinemetric fitting on 1-0 S(l) is limited both by lower S/N 
and smaller spatial extent of the data. The most notable 
differences between the velocity field of molecular gas and 
ionized medium is that the 1-0 S(l) is redshifted by 80 km 
s~^ relative to Paa at the nucleus and blueshifted by ~90 
km s~^ in the Northwest. Despite these differences, the H2 
velocity field is consistent with rotating disk with similar 
kinemetric parameters than derived for Paa (Fig. [7|. Rela- 
tively large Vobs — Vmodei residuals upto ^-^-100 km s~^ are 
detected Northwest of the nucleus coinciding with the re- 
gion where 1-0 S(l) is blueshifted relative to Paa (Fig. [8|). 
Interestingly, 1-0 S(3) velocity field is significantly different 
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Figure 8. Kinematic model of Pao; and 1-0 S(l) emission in the Super-antennae. The and /i4 components from Gauss-Hermite 
expansion of the Pao; profile are shown upper left and upper right, respectively. The difference between v^hs cii^d v^odel constructed 
without including higher order harmonics (rotation only) is shown for Pao; (lower left) and 1-0 S(l) (lower right). The contours in hs 
and h4 figures show the Pao; regions with Vobs - ^model = =t25 km s"-*^. 



from 1-0 S(l) northwest of the nucleus. The blue component 
of S(3) emission line is stronger than in S(l); thus the ex- 
citation temperature of the blue component is higher. This 
difference is unlikely to be due to extinction, as the wave- 
length difference between the lines is relatively small and 
therefore they both suffer from similar amounts of redden- 
ing. 



4 CONCLUSIONS 

We have presented the results of near-IR integral field spec- 
troscopy of the southern galaxy of the Super-antennae. The 
main results obtained from this project are: 

• Paa arises in a warped disc with position angle of 
~330°and inclination i=40-55°. The kinemetric parameters 
derived for H2 are similar. 



• Unlike in previous studies, we found the H-hand con- 
tinuum to be consistent with pure stellar emission, with no 
evidence for diluting non- stellar continuum at present. 

• The Paa/H2 line ratio maps reveal the presence of the 
ring- like structure with diameter ^-^1.5-2.0^^ coinciding the 
the region with lowest H2 ortho/para ratios. This morpho- 
logical argument together with non-thermal ortho-para ra- 
tios suggest UV fluorescence in dense clouds is an important 
H2 excitation mechanism in the nuclear region. 

• Extended H2 emission with r > lV3 has thermal or- 
tho/para ratio and very low Paa/H2 line ratios, indicating 
that thermal excitation in shocks is important at larger radii. 
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